Geometric Theory of Columnar Phases on Curved Substrates 
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We study thin self-assembled columns constrained to lie on a curved, rigid substrate. The curva- 
ture presents no local obstruction to equally spaced columns in contrast to curved crystals for which 
the crystalline bonds are frustrated. Instead, the vanishing compressional strain of the columns 
implies that their normals lie on geodesies which converge (diverge) in regions of positive (negative) 
Gaussian curvature, in analogy to the focussing of light rays by a lens. We show that the out of 
plane bending of the cylinders acts as an effective ordering field. 
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The concept of long range order is the lynchpin of 
much of our understanding of both phase transitions and 
the elasticity of the resulting crystalline and liquid crys- 
talline phases. In the presence of a quenched pertur- 
bation, this order can be disrupted plH- Moreover, in 
low dimensional systems thermal fluctuations can destroy 
long-range order: in two dimensions long range order can- 
not be achieved in systems with continuous symmetry 
Q, while in two and three dimensions one-dimensional 
periodic order is unstable to thermal fluctuations 
It would thus seem unlikely that there would be much 
to say about one-dimensional periodic order on a two- 
dimensional substrate. However, most of these "no-go" 
results suppose short range interactions. When the sub- 
strate is curved, however, long-range interactions aris- 
ing from global geometric constraints can lead to order- 
ing where none is usually expected. This feature opens 
up a novel venue for directing the self-assembly of thin 
layers of material with a desired structure by appropri- 
ately engineering the underlying substrate. Experimental 
realizations include self-assembled monolayers of tilted 
molecules Q, nematic double emulsions Q and solid 
rafts on curved interfaces 0]. The presence of geomet- 
ric constraints causes strong elastic deformations from 
the ground state of these materials that often cannot 
be captured by perturbative methods. In many cases, 
the investigation of these challenging issues of non-linear 
elasticity can benefit from blending classical continuum 
mechanics with modern geometrical methods @, |^, ■ 
The deformations of thin elastic shells is one such prob- 
lems that occupies a central place in any exposition of the 
basic principles of applied mechanics, but still continues 
to stimulate a wealth of theoretical and experimental re- 
search in analogous systems broadly defined as curved 
solids 0, 03. 

In this article, we address the less understood case of 
thin films composed of nearly incompressible columns 
such as those obtained from the self-assembly of block 
copolymers in a cylindrical phase fl3j . This system 
presents a natural arena to study how the curvature of 
the substrate affects the translational and orientational 
order that arises by minimizing the distortion energy of 
the cylinders [lij]. Cylindrical domains of block copoly- 



mers also have technological potential as nanowires and 
templates for quantum dots [15]. Non-planar confining 
geometries, for example cylinders with radius compara- 
ble but incommensurate to the copolymer period, have 
been successfully employed to generate structures not ac- 
cessible under planar confinement In the opposite 
limit, films of block copolymer cylinders (with tens of 
nm radius) have been synthesized on substrates made of 
larg er, melted colloidal particles (with few pro. radius) 
\u\ . There is markedly increased orientational order of 
the cylinders on these curved substrates when compared 
with flat substrates. Here we propose a coupling mech- 
anism between the local normal to the columns and the 
varying curvature of the substrate that acts as an effec- 
tive ordering field for the columnar order. In particular, 
we argue that where the Gaussian curvature of the sub- 
strate changes sign a boundary-like term, akin to a weak 
anchoring term in a liquid crystal cell, aligns the columns. 
The constraint of equal spacing then dictates the order 
in the region of negative Gaussian curvature. 

The free energy density for columns has two essential 
ingredients: the compression energy which sets the in- 
tercolumnar spacing, and the bending energy measured 
by the three dimensional curvature, n, of the columns. 
The curvature is further decomposed into the geodesic 
curvature, n g , and the normal curvature, n n : We denote 
the normal to the surface as N, the unit tangent vector 
of the columns as t, and define n = N x t to be the 
normal to the columns in the tangent plane of the sur- 
face. Then k 2 = \d s t\ 2 = |(t • V)t| 2 , n g = n • (t • V)t 
and K n = N • (t • V)t [18], from which it follows that 
The normal curvature measures the bend- 
ing of the cylinder as it is constrained to lie on the sur- 
face whereas the geodesic curvature measures the intrin- 
sic curvature of the cylinder bending within the surface. 
For example, great circles on a sphere have vanishing 
geodesic curvature yet nonzero normal curvature because 
they still bend. 

The compression energy density is most easily ex- 
pressed in terms of a phase field ^(x), which is related to 
the density p = po+Pi cos(27r$/a) of one of the polymeric 
blocks. Here, a is the preferred intercolumn spacing, po 
is the constant background density, and p\ is the order 
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circle of K=0 columns (solid) geodesies (dashed) 





































































L:-L, : square patch------ 




























1 


1 







-10 12 3 



FIG. 1: Configuration for columns (solid curves) shown in 
projection from above on a surface with height function 
h(x,y) = ho exp[—(x 2 + y 2 )/(2ro)]. We set the column at 
x = — lOOro to be straight and the red circle is the line with 
K = 0. The column normals lie on geodesies (dashed lines) 
that diverge and converge in regions of negative and positive 
Gaussian curvature, respectively. The dark lines delineate a 
triangle made of three geodesies: the x-axis, the straight col- 
umn at x = — oo, and the geodesic that bends towards the 
a;-axis. In A) we show the full bump for ho = 4ro. Note that 
some geodesies turn by 7r. In B) we depict the region y > 
for a bump with ho — 2ro. For clarity, we stop the geodesies 
when they intersect other geodesies. 

parameter for one-dimensional crystalline order. In this 
language, the columns sit at the level sets of <E> = ma 
(to G Z) and the columns are at their equilibrium spac- 
ing when | V$| = 1. The compressional strain, e, must be 
written in terms of | V$| and must vanish when | V<f>| = 1; 
for example, e = | [l — (V<&) 2 ] is one of the many ac- 
ceptable forms. The precise form of e is unimportant for 
our discussion, as we shall see. The total free energy in 
the columnar phase is 

F=| |^[e 2 + AV] , (1) 

where A = y/Kx/B is the penetration depth, K\ is the 
splay modulus, and B is the bulk modulus. In the 
experiments of Ref. [I?) A « 3 nm is a tenth of the 
cylinder spacing. It is useful to note that the deriva- 
tive of |V$| = n • V$ perpendicular to the columns 
is P/dj|V$| = |V$|(n • T))m where P'( = Sf - n^rii is 
the transverse projection operator and D is the covari- 
ant derivative in the surface. If the strain vanishes, then 
|V$| = 1 everywhere and (n • D)n = 0. We recognize 
this as the geodesic equation for n, and conclude that 
vanishing strain implies that the normals to the columns 
are geodesies on the surface. 

This observation suggests an analogy between this sys- 
tem and geometric optics, where the column normals take 
the place of light rays. The Gaussian curvature K of 



the surface plays the role of the index of refraction in 
focussing the geodesies, as is evident from the geodesic 
deviation equation, d 2 ^/ds 2 = —K£, where s is the ar- 
clength along the geodesies and £ is a measure (precisely, 
the Jacobi field which points from one geodesic to an- 
other [19() of the distance between two infinitesimally- 
close and parallel geodesies. Thus, K < causes the 
column normals to diverge while K > causes them 
to converge. We illustrate this in Fig. [1] for columns 
aligned by a smooth wall far to the left of a Gaussian- 
shaped bump, which sets up a family of initially parallel 
geodesies. The figure shows the effect of geometrical frus- 
tration on column-by-column growth that proceeds from 
left to right. The bending of the geodesies leads to sharp 
kinks in the columns downstream from the bump sim- 
ilar to the formation of caustics in ray optics. Indeed, 
this analogy between geometry and optics is the basis 
for gravitational lensing [2(| and the columnar system 
provides a "down to earth" realization of this important 
observational tool. Note also that geodesies can turn by 
7T and return to the column on the left if the bump has 
a large enough aspect ratio. 

The angle a which describes the downstream grain 
boundary where two families of columns with different 
orientations meet can be related to the Gaussian curva- 
ture of the surface by utilizing the Gauss-Bonnet theorem 
on a surface patch M [l9| : 

I K g ds = 2tt - V A6i - / KdA , (2) 

JdM i JM 

where A9i are the external angles at any vertices in dM , 
and the integral around dM is broken into piecewise in- 
tegrals along the smooth part of dM. For the geodesic 
triangle depicted in Fig. [TJ this calculation yields 

a = / KdA , (3) 

J M 

since the three edges are geodesies, two of the external 
angles are 7r/2 and the remaining external angle is n — a. 
We take the left edge of the triangle out to x — > — oo. 
Notice that as the distance from the center of the bump 
increases, the triangle covers a larger portion of the sur- 
face, asymptotically containing up to half of the bump. 
We conclude, therefore, that a approaches zero with in- 
creasing x since the integrated Gaussian curvature of any 
smooth bump tends to zero in this limit. We expect that 
near the kinks, the equal-spacing constraint will relax so 
as to lower the overall energetics of the caustic structure. 
Whether dislocation pairs unbind along the y-axis in Fig. 
[T] to soften the curvature or whether other low energy 
configurations Q which share the same topology come 
into play is an open question. In contrast to the case of 
crystals on curved substrates [ll|, [12], there is no local 
obstruction to growing equally spaced columns. Given 
any single column, the equal-spacing condition |V$| = 1 
can be integrated to determine the shape of the remain- 
ing columns within a finite size region. The compression 
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strain does not, therefore, impose any local orientational 
order on the columns. 

Since at every point on the surface there are geodesies 
in every direction, one might think that the columns 
could lie along geodesies. However, this assumption is 
incompatible with equal spacing on a curved substrate. 
To see this, assume that the columns lie along geodesies 
and that they are equally spaced. We start with one 
geodesic column. At two points on the column L t apart, 
we move off along the normal direction and trace the 
geodesies obeying (n • D)n = 0. After integrating both 
a fixed distance L n , these geodesies will end on a new 
column and n will still be perpendicular to the column's 
tangent vector. Here again, the Gauss-Bonnet theorem 
(Eq. [2|) constrains possible column configurations. In the 
case of our L t x L n rectangle (shown in Fig. [T]), there are 
four vertices with A9 = ir/2 and so 

f K g ds = - KdA . (4) 

J dM JM 

However, we assumed that n g — along the normal flows 
and along the columns. We see this is only possible for 
arbitrary rectangles if K = on the substrate! We have 
thus established the incompatibility of equal spacing and 
vanishing geodesic curvature of the columns. The true 
ground state will be some compromise between a non- 
ideal spacing and some curvature of the columns (the 
configuration shown in Fig. [I] depends on the growth 
history and is unlikely to be the actual ground state). 
However, we note that it is the gradient of the compres- 
sion that is proportional to the curvature of the normals 
and so, dimensionally, the compression energy will grow 
with an extra factor of L 2 (in a region of size L) com- 
pared to the bending energy. We thus expect, at long 
distance scales, that the curvature of the columns will 
absorb the effect of the Gaussian curvature and that the 
column spacing will be set at the ground state value, i.e. 
e = 0. 

Based on this discussion, it would be natural to sup- 
pose that, since the lines of longitude are geodesies on 
a rotationally symmetric bump, the columns lie along 
lines of latitude. However, the curvature of the columns 
frustrates this expectation. The intrinsic curvature does 
not give us a local way of selecting one configuration 
over another - one must compute the total curvature 
energy for different configurations with e = (Includ- 
ing defects in the ground state, along the lines of Ref. 
[l2| is beyond the scope of this paper). We turn, there- 
fore, to the extrinsic curvature of the columns, K n = 
N • (t • V)t = —titjdjNi, where the final equality follows 
from N ■ t = 0. If we express t = cos#ei 4- sin#e2 in 
terms of the principal directions e^, the extrinsic curva- 
ture becomes n n = K\ cos 2 9 + K2 sin 2 9 where K>i £1X6 the 
principal curvatures [l9j]. We take ei to be along the ra- 
dial direction and e2 along the azimuthal direction. The 
minimum of k 2 depends on whether K = k\K2 is positive 
or negative: when K < 0, the two principal curvatures 
differ in sign and k 2 attains its minimum when K n van- 




FIG. 2: A) Data from [T^]. Dark lines are the columns. The 
central region has K > 0, and the outer region has K < 
0. B) + C) Equally-spaced columns lying on geodesies of a 
cylindrical surface, a surface with K = 0. In B) the columns 
are straight in three-dimensions and have c„ = 0, while in 
C) the columns are curved in three dimensions and K n 7^ 0. 
Near the circle of K = 0, the bump shares the geometry of 
the cylinders. 

ishes, i.e. tan 2 #o = —Ki/k.2- Deviations 5 = 9 — 9q 
from the preferred angle have an energy per unit area of 
Af = -2B\ 2 K6 2 + 0{S 4 ). 

For simplicity, we consider a substrate which is a sur- 
face of revolution with height function h(r). Our argu- 
ments do not depend on the particular shape of the bump 
but rather on the generic fact that near the top K > 
and at the bottom, in the "skirt", K < 0. Since we want 
to consider smooth, azimuthally symmetric bumps, we 
know that h'(r) vanishes at r = and r = 00 and so 
somewhere in between, at r = ro, h"(ro) — 0. Since the 
radial direction is also a principal direction, the princi- 
pal curvature K r vanishes there, and subsequently K = 
on that circle. On this circle, we have a universal, pre- 
ferred direction for the columns - they lie along the ra- 
dial direction since tan 2 9q = 0. Thus the circle at radius 
r = ro imparts a boundary condition, albeit a very soft 
one. Near r — the restoring energy is quartic in 5, not 
quadratic. A radial alignment of the columns appears to 
be the energetically favored solution on the ring r = tq for 
sufficiently steep bumps. This picture is roughly borne 
out by experimental observations reproduced in Fig. [2] 
which shows a highly disordered columnar arrangement 
near the top of the bump, a radial corona of aligned 
columns near the ring r = ro, with a dramatic transi- 
tion to columnar rings at slightly larger radius. 

Focussing now on the top of a bump, where K > and 
r < ro, the minimum of k 2 occurs when sin20o = 0. If 
«2 > «i then the minimum is at #0 = with a deviation 
energy of Af ~ ^\ 2 {n2 — ki)ki5 2 , while when K2 < Ki, 
the preferred angle is 9q — tt/2 with a deviation energy of 
Af ~ -jA 2 (ki — K2)k2<5 2 - Since n\ grows from zero near 
ro, the radial ordering will still be favored near the band 
of vanishing K. However, depending on the details of the 
bump, there will either be azimuthal or radial ordering 
near the top. The incompatibility of radial ordering in 
the K > region would force dislocations and disclina- 
tions into the texture since radial aligned columns have 
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Kg = and so, by our original arguments, equal spacing 
would be impossible. It may be difficult, however, to dis- 
tinguish between a slightly disordered azimuthal ordering 
(«i > K2) and a dislocation texture arising from bound- 
ary conditions at r — tq when K2 > K\. In contrast, 
the columns on the flanks of the bump are normal to ra- 
dial geodesies as we may have expected from our initial 
discussion. This extrinsic effect should be contrasted to 
the intrinsic coupling between the column direction and 
gradients of intrinsic quantities considered in 12 If . 

As stated above, the data of Hcxemer et al. [17| finds 
azimuthal ordering in the K < regions of the bump 
with less pronounced order in the K > regions. More- 
over, in that experiment there is no long range order in 
the flat portions of the sample where the restoring force 
proportional to K vanishes. The translational correlation 
length between cylinders does not exceed 200nm, corre- 
sponding to a patch of less than ten equidistant columns. 
On larger distances, the system appears to be comprised 
of orientationally uncorrelated patches even if a configu- 
ration of parallel cylinders is the obvious ground state for 
perfectly annealed flat samples. These results are consis- 
tent with prior experiments on flat substrates [22| . 

Though the asymptotic shape of the bump controls 
the asymptotic lines of vanishing normal curvature, the 
rigidity of equal spacing in the K < region and the inde- 
fatigable radial ordering close to the inflection point pro- 
vides (in the absence of defects) a mechanism that gen- 
erates involutes along the surface . This follows from 
the discussion of the geodesic curvature: the violation of 
equal spacing is more energetically costly than deviations 
from the preferred curvature. In the experiments [l7j , an 
array of bumps were used as a substrate instead of an 
isolated bump, as we have considered here. In this situ- 
ation the hexagonal geometry implicitly generates more 



boundaries which act as ordering fields which are likely to 
promote columnar ordering in the skirts. The texture of 
bumps acts as a source for large core disclinations in the 
normal field [l(| • Our analysis considers a single layer of 
columns where the bending along their normal direction 
was neglected. The experiments imaged the top layer of 
a film roughly ten layers thick which was grown on the 
curved substrate. Unimaged dislocations and strains are 
necessary below the surface to accommodate the defor- 
mation of the columnar crystal and are formed as the 
layer grows on the curved surface. How the underlying 
layers couple to the two-dimensional geometry of the top- 
most columns is an open question. 

In summary, we have studied how the curvature of a 
substrate influences the order of incompressible, periodic 
structures by drawing an analogy to geometrical optics. 
Regions of positive Gaussian curvature cause the con- 
vergence of the column normals, negative Gaussian cur- 
vature lead to their divergence. On the one hand, the 
extrinsic curvature of the columns results in a local, pre- 
ferred ordering depending on the shape of the surface. 
On the other hand, on a symmetric bump the region near 
K = provides a generic mechanism of radial alignment 
akin to a soft boundary condition. We hope to discuss 
the role of the breakdown of the equal spacing condition 
and the resulting dislocations in future work. 
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